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Abstract 

We discuss the power scaling rules along the lines of a complete Heavy Quark Effective Field 
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1 Introduction 



The application of effective theories directly derived from first principles - Quantum 
Chromodynamics (QCD) in particular - to heavy flavor physics has currently become 
a very important tool to cope with the complexity of the strong interaction dynamics. 
One of its advantages is the appearance of symmetries in the very large quark mass 
limit, not manifest in the full theory. Moreover, symmetry breaking effects due to the 
finiteness of quark masses can be systematically incorporated via power expansions 
in the reciprocal of the heavy quark mass. One of the foremost applications of 
heavy-quark spin-flavor symmetry has led to accurate and almost model-independent 
extractions of some Kobayashi-Maskawa matrix elements in inclusive and exclusive 
B meson decays (see for a review). 

Over the nineties, a heavy quark effective theory (HQET) was developed |2|, |3| 
to deal with the phenomenology of heavy-light mesons and baryons. Moreover, 
there were early attempts to extend its initial scope to heavy-heavy bound states 
like the J/ip resonance |4| or the B c meson and doubly heavy baryons ||. On 
the other hand, a close approach also emerged in parallel, especially designed to 
describe heavy quarkonia inclusive decay and production, the Non-Relativistic QCD 
(NRQCD) || [7|. In this formalism, the decay widths or cross sections are factorized 
as a sum of products of short- distance coefficients and long-distance NRQCD matrix 
elements. 

Although both HQET and NRQCD Lagrangians present a formal analogy, there 
is a basic difference between both frameworks. In the former one, there are basically 
two relevant scales: the heavy quark mass tuq and Aqcd, the characteristic parame- 
ter of the strong interaction. In NRQCD, besides Aqcd, there are (at least) another 
two low-energy scales: rriQv' (soft) and rriQv' 2 (ultrasoft) []. 

Recently, several applications of HQET involving both heavy quark and antiquark 
fields altogether to the study of the phenomenology of the weak and strong interac- 
tions have appeared in the literature (|8|, [| and [O, respectively). At the same time 
it becomes suitable a careful survey of the foundations on which such a theoretical 
framework is based. Therefore, in this paper firstly we shall focus on the quantum 
aspects of the heavy quark effective fields, summarizing some work presented in pre- 
vious references [[□], [EJ. In particular, we shall keep both particle-antiparticle mixed 
sectors of the effective Lagrangian, in accordance with the pioneering work by Wu 



in Ref. [|13|], further developed and discussed in ||14j| . Hereafter we shall refer to this 



extended framework as HQEFT as in [|J, to distinguish it from "usual" HQET. 

Furthermore, one of the goals of this paper is to analyze whether HQEFT could 
provide an appropriate framework to describe meaningfully the hadropro duct ion of 
charmonium. In particular, we shall analyze the power counting of four-fermion local 
operators, whose scaling rules for charmonia have been recently hypothesized in Ref. 



15] along the lines of HQET. We shall examine this issue with the help of diverse 



extractions of some matrix elements performed by several authors. 



1 In this paper we shall employ v' to denote the three- velocity of heavy quarks inside quarkonium while v will stand 
for the hadron four-velocity, as usual in HQET; since residual momenta of heavy quarks will not be neglected in the 
hadron rest frame, v' obviously does not exactly coincide with the spatial component of v, i.e. v' |vj. 
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2 Notation and definitions 



Let us start following Ref. [11] by writing the plane wave Fourier expansion for a 
fermionic field Q v (x) corresponding to an almost on-shell heavy quark or antiquark 
(we omit indices relative to flavor and color) moving inside a hadron with four- 
velocity v, 

Qv{x) = Ql + \x) + Q[-\x) = I & £ [ &r ( p ) Ur ( p ) e -4« + fot(p) ^(p) e *« ] 

r 

(1) 

where r refers to the spin and J stands for the chosen normalization; 6 r (p)/6j(p) is 
the annihilation/creation operator for a heavy quark/ antiquark with three- momentum 
p (p° ~ + y/m 2 + p 2 ) . Let us firstly focus on the particle sector of the theory. 

Particle Sector 

As the heavy quark is almost on-shell one should require for each Fourier com- 
ponent that p 2 = trig + A 2 where A 2 satisfies lim^^co A 2 /ttvq — > 0. From now on, 
we shall set A 2 = as a good approximation for heavy quarks bound in hadrons. 

On the other hand, spinors are normalized such that ul(p)u s (p) = 2p°N 5 rs and 
the creation/annihilation operators satisfy: 

[b r (p'),bl(p)]+ = K5 rs 5 3 (p , -p) (2) 

where K, N are the corresponding normalization factors. 

Now let us redefine the momentum of each Fourier component in Eq. (1) according 
to HQET as the sum of a mechanical part and a Fourier residual four-momentum k, 

p = ijiqv + k (3) 

Hence one can write from (1) 

— J2 b r (v,k) u r (k) e~ ik -* (4) 

r 

The main point to be stressed again is that we shall require that each Fourier 
component should satisfy the almost on-shell condition Pq ~ nig. Therefore, 

k 2 + 2m Q vk = (5) 

Thus, in the hadron rest frame, the kinetic energy of the heavy quark v-k = k° is 
related to k through the constraint 

(k ) 2 + 2m Q k° - k 2 = (6) 



which yields the expected relation for the positive root k° = \/m>Q + k 2 — rriQ. In 
the non-relativistic limit obviously, k° ~ k 2 /2m,g. 
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Notice that the annihilation (and creation) operators and spinors have been simply 
relabeled in Eq. (4), satisfying the same normalization as above, though expressed in 
terms of the Fourier residual momentum k. In particular, b r (v, \z)/b\(v, k) correspond 
to annihilation/creation operators of a heavy quark with residual momentum k in a 
hadron moving with four-velocity v, satisfying accordingly 

[b r (v, k'), bl(v, k)] + = K 6 rs 5 3 (k' - k) (7) 

and 

Y, r u r (k)u r (k) = N [ m Q (l + f) + # ] (8) 
where the normalization factors must obey the combined relation: [T(| ^ 



KN 1 1 

J 2 ~ (2tt) 3 2p° ~ (2vr) 3 2( y m Q v° + k°) ^ 

On the other hand, it is quite usual in the literature to identify effective heavy 
quark fields with those "leading" components of the Fourier expansion corresponding 
to momenta p M ~ vtlqv^ (or equivalently with fc M components close to zero), so they 
rather look like single spinors or anti-spinors at leading order in l/mg. However, 
this is a too restrictive view since in constructing HQEFT one should allow the k 
components of the effective quantum fields to range over values of the order of Aqqd- 



Next let us introduce the effective fields in the following standard manner || [17 



h (+)( x ) = e im *™ l±i Qi + \x) (10) 

HW(x) = e im ^ x ^ Qi + \x) (11) 

where Hv (x) represents the "small" (lower in the hadron reference frame) compo- 
nent of the corresponding spinor field. Thus one may identify from the expansion 
(4) 

h£\x) = ^ / E b ^ u r(k)e-^ (12) 

r 

Hi + \x) = ^ / ^ E b r (v,k) u r (k) e~ ik -* (13) 

r 

Sometimes it can be read in the literature that annihilates a heavy quark 
whereas allegedly Hi + ^ creates a heavy antiquark, both with velocity v. Although 
may take into account the antiparticle features of a relativistic fermion, its 
Fourier decomposition does not allow for such a process, as can be seen from (13). 



2 Note that /^l/(27r) 3 2/fc° since k and p are not related through a Lorentz transformation; indeed, Eq. (3) represents 
an energy-momentum shift, to get rid of the unwanted frequency modes to be subsequently eliminated in the effective 
theory. 
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3 Anti-particle Sector 

Proceeding in a parallel way as in the particle sector, let us remark however that now 
the Fourier expansion of Qi~\x) involves negative frequencies. Therefore, starting 
from Eq. (1) we introduce the Fourier residual momentum in this case as 

p = m Q v - k (14) 

so k° will explicitly exhibit its negative sign. In effect, let us write 

— ~ b K k ) V r( k ) e ~ tk ' X ( 15 ) 

r 

The slightly off-shellness condition p 2 ~ rriq now implies from Eq. (14): 

k 2 - 2m Q vk = (16) 
which can be written in the frame moving with velocity v as 

(A; ) 2 - 2m Q k° - k 2 = (17) 

whose negative root reads k° = rriQ — ttIq + k 2 ~ — k 2 /2mg. 
The effective fields are then defined as 

hi~\x) = e -*»Q- izi Q(-)( X ) = ^ / ^ E « r (k) e-**"* (18) 

r 

HtHx) = e-^ ±±t Qi'Hx) = l -±t [ ^ J2 blMvrMe-*" (19) 

r 

where the anti-spinors satisfy 

J2 r v r (k)v r (k) = N [ m Q {j - 1) - # ] (20) 

4 A complete tree-level HQEFT Lagrangian 

Once introduced the notation and definitions for the heavy-quark effective fields in 
the previous Section, we actually get started by expressing the Lagrangian in terms 
of the effective fields keeping all non-null terms, leading in principle to the possibility 
of describing annihilation or creation of heavy quark- ant iquark pairs at tree level. 
Indeed after performing an energy-momentum shift by introducing a center-of-mass 
residual momentum, only low energy modes of the fields (about the heavy quark 
mass) would be involved making meaningful our approach within the framework of 
an effective theory. 
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Thereby a crucial difference of this approach w.r.t. other standard works is that 
we are concerned with the existence of two-fermion terms in the transformed La- 
grangian mixing large components of the heavy quark and heavy antiquark fields, 
i.e. hl^Th^, where T stands for a combination of Dirac gamma matrices and covari- 
ant derivatives. The implications on the description of heavy quarkonia production 
will be seen in Section 6. 

The tree-level QCD Lagrangian is our point of departure: 



£ = Q (ip - m Q ) Q 



where 



Q = Q(+)+Q(-) = e -im Q vx h (+) +H (+) 



_|_ e im Q v-x 



(21) 



(22) 



and D standing for the covariant derivative 

^ = _ igTaA » 

with T a the generators of SU(3) C . Substituting (22) into (21) one easily arrives at 

C = £(++> + + £("+) + £(+") (23) 

where we have explicitly splitted the Lagrangian into four different pieces correspond- 
ing to the particle-particle, antiparticle-antiparticle and both particle-antiparticle 
sectors. The former one has the form 

£ (++) = h(+) iv .^h^ - H ( f\iv~L) + 2m Q )H^ + h^0 ± H^ + H^0 ± h^ 

(24) 

corresponding to the usual HQET Lagrangian still containing both hi + ^ and 
fields. We employ the common notation where perpendicular indices are implied 
according to 

D± = D a (g^-v»v a ) 
Regarding the antiquark sector of the theory 

= -U- ) iv~l}h { - ) + H(-\ivTt - 2m Q )H i f ) + h^ifx.H^ + H^ifji^ 

(25) 

The latter expressions, considered as quantum field Lagrangians, certainly do not 
afford tree-level heavy quark-antiquark pair creation or annihilation processes stem- 
ming from the terms mixing hl ±] and fields since they contain either annihila- 
tion and creation operators of heavy quarks or annihilation and creation operators 
of heavy antiquarks separately. 

Nevertheless there are two extra pieces in the Lagrangian (23): 

£(-+) = e -i2m Q vx x ( 26 ) 
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and 



£(+-) = e i2m Q vx x ^ 



where use was made of the orthogonality of the and Hv fields flifl . As could 



be expected, there are indeed pieces mixing both quark and antiquark fields leading 
to the possibility of annihilation/creation processes. After all, Lagrangian (23) is 
still equivalent to full (tree-level) QCD []. 

Let us remark that, at first sight, one might think that the rapidly oscillating 
exponentials in Eqs. (26-27) would make both £(~ + ^ and C^ + ^ pieces to vanish, 
once integrated over all velocities according to the most general Lagrangian |T8" 



However, notice that actually this should not be the case for momenta of order 
2rriQV of the gluonic field present in the covariant derivative. In fact, only such 
high-energy modes would survive, corresponding to the physical situation on which 
we are focusing, i.e. heavy quark-antiquark pair annihilation or creation. 

From Eq. (26), the heavy quark-gluon coupling for an annihilation process is 
given by 

r(~+) _ ~i2m Q v-x rp A a 

'-coupling ~ e y^a^V * 

[ (a) 

- H ( v +) (6) 

+ hi- } ^h^ (c) 

+ H^^H^] (d) (28) 

There is an equivalent expression coming from jC^ u J ling just by means of the 
substitution v— > — v, corresponding to a heavy quark-antiquark pair creation process. 

Let us make an important remark concerning the gluonic field. As can be seen 
at once, the exponential factor e ~ l2m Q v - x i n Eq. (28) cancels against the strong x 
dependence of the A a ^ field creating a hard gluon. Therefore we can write 

A s ° ft (x) = e~ i2m Q v - x A^x) (29) 

Therefore, once removed the strong dependence of the gluon momentum on 
(twice) the heavy quark mass, we are left with an effective gluonic field A 8 "'*, with 
a mild dependence on the residual momentum, similarly as for heavy quark effective 
fields. Thus, the physical description of the annihilation process shown in Fig. 1 
could be performed in terms of massless quarks and a soft gluonic field. This is pretty 
similar to the situation reached in heavy-light systems, constituting a cornerstone in 
our approach to be developed in the following Sections, providing a theoretical basis 
for the use of HQEFT to describe heavy quarkonium production. 



Superscript "(+)/(—)" on the effective fields labels the particle/antiparticle (i.e. heavy quark/antiquark) sector of 
the theory Actually hi + ^ (hi~^) corresponds to negative (positive) frequencies associated with creation (annihilation) 
operators of quarks (antiquarks). In fact some extra " + /— " signs should be added on the conjugate fields, i.e. 
and Hi~^ + , which however will be omitted to shorten the notation. 
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2mv+q-q' 



Figure 1: Heavy quark-antiquark pair annihilation into a gluon of momentum Ittiqv + q — q' . In the center-of 
mass system q = (q , q) and q' = (—qo, q). This process can be described meaningfully by HQEFT if q 2 << tUq, 
i.e. the square invariant mass of the gluon is close to 4mg. 

Next we want to eliminate the unwanted degrees of freedom associated to the 
"small" components Hv in (a), (b) and (d) in Eq. (28). Notice that the piece 
labeled as (c) is the leading one in the above development. 



5 Derivation of the annihilation vertex 



In all our later development we shall assume almost free quarks in the initial (or 
final) state (see Fig. 1). Hence, heavy quark fields appearing on the external legs 
of the Feynman diagram should be solutions of the unperturbed Dirac equation of 
motion. Hence we can write 



V 
V 

and for the conjugate fields 

V 

BP - 



(iv- d + 2m Q - ie)~H $ ± h { v +) 
(—iv-Jf + 2vrtQ — (j) 



hl + \ (j) j_(iv -*W — IrrtQ + it) 1 
h^i (}) ±(—iv-*~d — 2itlq + ie)^ 1 



(30) 
(31) 

(32) 
(33) 



Let us note that we are using the free particle equations which can be viewed as 
derived from the non-interacting parts of the Lagrangians C^ ++ ^ and & -* respec- 
tively, i.e. we neglect the soft gluon interaction among heavy quarks amounting to 
describe them as plane waves, i.e. actually no bound states as a first approximation 



In momentum space we will write 
u(m Q v + q) = 



1 + 



2mg + v-q 



Uh(q) 



(34) 



v[m Q v - q 



Vh{-q 



1 - 



2m Q - v-q' 



(35) 



where u(p) denotes the full QCD spinor whereas Uh(q) represents the HQET spinor 
i.e. obeying jju^ = Uh\ similarly v(p) denotes the full QCD antispinor and Vh(—q')) 
represents the HQET antispinor obeying jjVh = —Vh- 
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Actually, in order to arrive to Eqs. (34-35) from Eqs. (30-33) denominators have 
to be expanded as power series of derivatives acting on the x-dependent factors, 
assumed exponentials, to be finally resummed as geometric series of ratio v-q/2rriQ = 
—v-q'/2mQ = —q 2 /4:mQ according to the on-shell conditions (5) and (16). As a 
consequence, Eqs. (34-35) are only valid under the condition — q 2 < 4mg, which 
implies q 2 < 8rriQ. Therefore, the requirement q 2 < niq satisfies the above condition, 
allowing a non-relativistic expansion. 

From the &~ + ' Lagrangian we readily get for the on-shell heavy quark (vector) 
current coupling to a gluon, suppressing color indices and matrices, 



+ 



-v fl + 



2m Q + v-q (2m Q + v-q) 2 _ 



Uh 



(36) 



which also can be written as 



11 H ~ : 1- 



2m Q + v-q (2m Q + v-qy 



(37) 



since 



14 



with the projectors P± = (l±^)/2, and q' ±u v u = qj_ u v u = 0. 

Equations (36-37) are Lorentz covariant and hence valid in any reference frame. 
In the following Section we focus on the center-of-mass system, where some simpli- 
fications occur. 



5.1 Center-of-mass frame 

We shall make use of the anticommutation relation: 

Now, in the center-of-mass frame we can write 
pression for the heavy quark vector current 

2E, 



v h 



7± + 



E q + m Q ' x (E q + m Q ) 2 _ 



Uh 



-q 2 , leading to the ex- 



(3f 



where E q = uiq + q° = yrrig + q 2 . Moreover, the following identity is satisfied (in 
the Dirac representation): 



(39) 



Therefore identifying Uh = ^jE q + rng £, and Vh = ^jE q + mq r], we obtain from 
(38-39), 



2E a o 3 + 



2g J cr-q 

{E q + m Q )_ 



(40) 
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i.e. the vertex obtained from full QCD in terms of the Pauli spinors £ and 77, allowing 
a systematic non-relativistic expansion: the leading term 2m(j]'a^) associated with 
(c) in Eq. (28) as expected. 

Had we started from the Lagrangian i.e. focusing on a creation process, 

the heavy quark vector current would be 



2E„ 



and since 



we rapidly recover 



E q + m Q 

p+ L p- 



2E q a 3 



+ 



2q 



"-< . 



(Eq + m Q f 



Vh 







2q J cr-q 



(41) 



(42) 



(43) 



(E q + m Q ) 

which is the analogous expression to Eq. (A. lib) in the c.o.m. frame as given in 



20 1 . 



In sum, we have derived in this Section the vector current of on-shell heavy quarks 
coupling to a hard gluonic field directly from the corresponding Lagrangian h \ 
pieces, giving consistency to our formalism. 



6 Heavy Quarkonia Hadroproduction 

In 1992-93, CDF and DO Collaborations surprisingly found an excess of charmonia 
prompt hadroproduction at the pp Tevatron collider |2l[ w.r.t. to the theoretical 
expectations based on what was considered at that time as conventional wisdom, 
the so-called color-singlet model (CSM) pifl . This disagreement was particularly 
amazing because data spread over large transverse momenta where the theoretical 
analysis ought to be quite clean. Therefore another mechanism, known as the color- 
octet model (COM), was proposed by Braaten and Fleming in 1995 ]2"3|| , superseding 
(generalizing actually) the CSM. Soon later, the COM was viewed as deriving from 
NRQCD 0, giving an adequate framework for the factorization of the inclusive 
production cross section into short- and long-distance parts, and providing a sound 
basis to the theoretical analysis so far done. 

However, current problems in the interpretation of charmonia hadroproduction 
30| , especially regarding the (expected but unobserved) transverse resonance polar- 



ization, have cast serious doubts on the validity of NRQCD when applied to these 
processes. Actually, the meaningful application of the effective theory relies on the 
presumed convergence of the expansion governed by the typical velocity of heavy 
quarks v' and the strong coupling constant a s . Different contributions are assumed 
as leading whereas others are neglected under the assumption of the validity of the 
velocity scaling rules derived from NRQCD, mainly based on dimensional counting 
and very general physical arguments. 
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In this context some authors |T5| have recently hypothesized that the correct 
power counting for charmonia should be the dimensionless parameter KqcdI^c-, 
along the lines of HQET [], leading to predictions which could differ from the ex- 
pectations coming from the usual velocity scaling rules of NRQCD. In the following 
and based on our previous development, we want to give a theoretical support to 
such a point of view. In fact we are advocating in this paper that HQEFT could 
provide an appropriate and consistent framework for the analysis of hadropro duct ion 
of heavy quarkonia. To this end, from Fig. 1 we may identify the heavy quark three- 
momentum q relative to the c.o.m. as the residual three-momentum k introduced 
in Section 2. 

On the other hand, in charmonium systems it seems sensible to assume v' ~ \q\/rriQ 
~ A/rriQ, where A stands for a "typical" scale characterizing the quarkonium soft dy- 
namics, of order few times Aqcd- Therefore dynamical (i.e. non-Coulombic) gluons 
in cc bound states should be of the type (A, A) though the typical residual four- 



momenta of bound heavy quarks should be of the type (A /2m,Q, A) (see |24|, [31 
in close analogy with heavy-light hadrons. Hence the heavy quark field expansions 
given in Eqs. (12) and (18) about the heavy quark mass, and the reinterpretation 
Q f e -i2m Q v-x j^, j n ^29) as a soft gluonic field coupling to a heavy quark vector 
current make sense altogether. 

On the other hand, the effective QCD Lagrangian || to deal with heavy quarko- 
nium is usually written as: 

C eff = C 2 + A + Clight (44) 

where C 2 contains bilinear operators, £4 stands for the four-fermion piece and Cu g ht 
takes into account light quarks and gluons. 

In NRQCD heavy quarks are treated according to a Schrodinger-type field theory 
with separate two-component Pauli spinors (commonly denoted as ip an d x)i rather 
than with four-component Dirac fields as in HQEFT. Relativistic effects of full QCD 
are reproduced through correcting terms bilinear in the quark field or in the antiquark 
field. However, mixed two-fermion operators involving both quark and antiquark 
fields in £ 2 corresponding to creation/annihilation of QQ pairs are excluded. 

On the other hand, in order to account for heavy quarkonium annihilation or 
creation, a set of four-fermion operators are thereby included in the piece £4 of 
the effective Lagrangian (44), representing local interactions whose short-distance 
coefficients have to be determined by means of a matching procedure to QCD as 



described, for instance, in |p0 | 



Actually, in this work we are advocating that the £2 piece of the effective La- 
grangian can be extended as 

C 2 = 4 ++) + 4~ _) + 4~ +) + 4 + ~ ] (45) 

where the £ 2 + ^ and C 2 + ^ terms can be identified with those coming from Eqs. (26- 
27), therefore taking into account QQ pair formation and annihilation at tree-level. 



4 The authors of [|L5| introduce the acronym NRQCD C which we identify with HQEFT applied to charmonium. 
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Figure 2: Graph representing a qq annihilation followed by the creation of a QQ[ 3 Si] pair, evolving into a final 
spin-triplet quarkonium color-singlet state by emission of two soft gluons via a double chromoelectric dipolar 
transition. The two dotted lines delimit the gluon-quark vertex to be described along the lines of HQEFT 
following our suggestion of including the description of the QQ creation within the realm of HQEFT. Moreover, 
local four-fermion terms in the Lagrangian accounting for heavy quarkonium production can be thought as 
derived from the ) and h ' pieces of the same HQEFT Lagrangian. 

Notice that, as far as the heavy quark mass dependence has been removed from the 
effective fields, only low-energy modes (of the order of A around rriQ,) remain as 
stressed before. Hence the situation turns out to be very similar to the standard 
application of HQEFT to heavy-light systems. On the other hand, we keep the £4 
piece in the effective Lagrangian, although written in terms of HQEFT fields. 

For defmiteness, let us consider quarkonium hadroproduction through the par- 
tonic channel: 



where H denotes a particular quarkonium state. (The diagrammatic representation 
for this production process yielding a vector resonance is shown in Fig. 2.) 

6.1 Color-octet matrix elements from HQEFT 

Under the assumption of factorization, the cross section can be written for inclusive 
resonance production as 



where da stands for the production cross section of a QQ pair in a definite color and 
angular-momentum state labeled by n. The short-distance interaction, represented 
by da in (46), can be accounted for by perturbative QCD, whereas the long-distance 
process is encoded in the vacuum expectation values of four-fermion local operators, 
i.e. < > 0. According to NRQCD, the relative importance of the various terms 
in the factorization is determined by the order in a s and some dimensionless ratios of 
kinematic variables in the short-distance factors, and by a hierarchy based on the or- 
der in v' of the matrix elements. Hence, the velocity scaling of the MEs is determined 



qq _ g * _> (QQ) 8 [ 3 S 1 ) - H X 




(46) 
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by the number of derivatives in the respective operators and the number of electric 
and magnetic dipole transitions needed to reach the final physical quarkonium state 
from the initial QQ pair produced at the short-distance process. 

Pictorially, the long-distance evolution lies on the right of the rightmost dotted 
line in Fig. 2, belonging to the realm of the effective Lagrangian £4 and providing 
the probability for a heavy quark-antiquark pair to evolve into a (definite) final 
quarkonium state. On the other hand, the subprocess g* — > QQ can be ascribed 
to the piece of the effective Lagrangian from our viewpoint. Therefore we 

shall use the hi + ^ and hi ^ HQEFT fields when rewriting the relevant operators. For 
example, the 8 ( 3 Si) 6-dimensional operator || can be expressed in the form of a 
current x current term (see figure 3), i.e. 

8 ( 3 5x) = fcH 7^ T a fcW (a^a H ) ft« 7±fi T a fcH (47) 

where a^(a^) stands for annihilation (creation) operators of heavy quarkonium H, 
such that 



a} H a H = ^\H + X><H + X\ 



x 



The above 8 ( 3 Si) operator can be seen as deriving from the corresponding piece 
of the £4 Lagrangian || now expressed in terms of HQEFT fields, 



= hi + hl T a W hi-h±v T a (49) 

m Q 

The physical interpretation of (47) is that of a local four-fermion operator which 
creates a slightly off-shell QQ pair in a 3 Sf^ state and destroys it again, analogously 
as in NRQCD. Its vacuum expectation value < 0|O 8 ( 3 Si)|0 > = < 8 ( 3 Si) > 
corresponds to the usual color-octet matrix element of NRQCD. Notice, however, 
that the current x current structure can be thought now as deriving from the £3 + ^ 
and ^ Lagrangian pieces, involving only low-frequency modes and providing 
consistency to all the theoretical framework. High-energy modes are encoded in the 
short-distance coefficients (F 8 ( 3 5'i) in this particular case), as usual in low-energy 
effective theories. (Nevertheless, the matching procedure, although formally the 
same as in NRQCD, might be seen as "built in" because the effective fields in the 
four-quark local operators follow from £2-) 

Hence the power counting of the effective theory should follow the lines of HQEFT, 
as far as rrigv' ~ A, and the expansion parameter can be identified as A/tuq. How- 
ever, let us remark that the intermediate colored bound state (QQ)$ evolving into 
final charmonium actually is a kind of hybrid state [25], whose "size" (and even na- 



ture) is not accurately known. Therefore, one could expect A to be of order few 
times Aqcd, namely 700-800 MeV, and typically (A/tjiq) 2 ~ 1/4, but with large 
uncertainties. 
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Figure 3: Since the formation of the quarkonium state H takes place over distances of the order of 1/A, the 
use of HQEFT is justified. Local interactions give rise to the Os( 3 'S'i) operator among others, in terms of 
fields coming from the Li piece of the Lagrangian. 



6.2 Power scaling rules 

Under the assumption that HQEFT could be a suitable effective theory to deal with 
charmonia production, we will follow a reasoning analogous to [U| ^4] , arguing that 
there could be a hierarchy of the matrix elements, different from "standard" NRQCD. 
Therefore we expect that a single chromoelectric dipole (El) transition should scale 
as A/mg at the amplitude level. Since two El transitions are required to reach a 

3 5'{ 1 ' color-singlet state from a 3 S[ 8 ' color-octet one, the resulting suppression factor is 
(A/rriQ) 2 . On the other hand, a single chromomagnetic transition (Ml) should scale 
as A/m,Q, because of the time derivative on the gluonic field strength-tensor, picking 
up a typical energy of order A. Thus, at the cross section level, the power countings 
are expected to be of order (A/m<g) 4 and (A/??7.q) 2 for a double El transition and a 
single Ml transition, respectively. 

Now, defining as usual the linear combination of MEs f] (where r is an adjustable 
parameter varying between 3 and 3.5 for charmonia hadroproduction) as 

Mf(<*> = < Oi {nS) CS ) > + r ^ (50) 

m Q 

the former one, i.e. < QSq) >, should become leading in the overall power 

counting of Mf^ ■ 



Let us now define the ratio 



M? (nS) 



R ^ns) = r (n = l,2) (51) 

<ot (nS \*s 1 )> 

Since HQEFT predicts a typical scaling (relative to the color-singlet contribution) 
of (A/m Q ) 4 for < Of n5) ( 3 ^) >, and (A/m Q ) 2 for Mf (nS \ one should expect 

RHnB) „ f^Y (52) 



5 The little difference in shape between the and 3 Pj 8 ' contributions as a function of transverse momentum does 
not permit independent fits and only a linear combination of them can be extracted from experimental data J2S^, ^] . 
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which could also be expressed as 1/V in terms of the more conventional power 
counting based on the typical velocity of heavy quarks, under the assumption that 

v' ~ A/l7lQ. 

Conversely, in standard NRQCD both matrix elements < 0^ nS \ 3 Si) > and 
Mf^ nS \ should scale similarly (as v' 4 w.r.t. to the leading color-singlet component). 
Therefore no hierarchy should appear, and one would expect R^( nS ) ~ 1. 

In the following we carry out a check based on the available experimental in- 
formation [] extracted by several authors from fits to Tevatron data on charmonia 
hadroproduction. In tables 1 and 2 we show the values of the color-octet matrix 
elements < Og ^ ( 3 Si) > and Mf^ for J/if> and if>' respectively, and the corre- 
sponding ratios R^^ nS \ 

Table 1: Values (in units of 1CT 3 GeV 3 ) of < Og ( 3 #i) > and M^ 15 ^ matrix elements (r varies 
between 3 and 3.5) and their ratios R^^ 1S \ obtained from Tevatron data on prompt J /if) inclusive 
hadroproduction. Error bars are only statistical. 



ME: < Ot^^SJ > 








m 


11.9±1.4 


45.4±11.1 


3.8±1.0 


m 


10.6±1.4 


43.8±11.5 


4.1±1.2 






3.3±0.5 


14.4±2.8 


4.4±1.1 




m 


6.6±2.1 


66±5 


10.0±3.3 






3.9±0.7 


66±7 


16.9±3.5 



Table 2: Values (in units of 1CP 3 GeV 3 ) of < O g ( 3 <Si) > and matrix elements (r varies 

between 3 and 3.5) and their ratios R^( 2S \ obtained from Tevatron data on prompt ip' inclusive 
hadroproduction. Error bars are only statistical. 



ME: 


<ot^ ) ( 3 s 1 )> 


Mf 2 ^ 










3.7±0.9 


7.8±3.6 


2.1±1.1 






1.4±0.3 


3.3±0.9 


2.4±0.8 


m 




5.0±0.6 


18.9±4.6 


3.8±1.0 


IN 




4.6±1.0 


17.7±5.7 


3.8±1.5 








4.4±0.8 


18.0±5.6 


4.1±1.5 



We realize that indeed a hierarchy seems to exist for < Of ln6, ( 3 S' 1 ) > and Mf (nb) . 
Although R^ 1S ^ varies in a somewhat wide range, a value equal or greater than four 
seems to be favored. On the other hand, the ratio R^( 2S '> turns out to be slightly but 
systematically smaller than R^^ 1S \ but close to four as well, i.e. compatible with 
the expectation v' 2 ~ A 2 /wig ~ 1/4, already mentioned at the end of Section 6.1. 

6 We limit ourselves to those references where MEs extractions of both J/ip and ip' are given; notice that distinct 
theoretical inputs may have been used, as for example different parton distribution functions. 
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In fact one should bear in mind that tR^^ nS ^ is not a fixed quantity for the whole 
charmonium family (nor the velocity v'). This remark is in agreement with the 
expected larger value of the typical velocity v' for higher n states in the charmonium 
sector or, equivalently, a smaller ratio tuq/A. 

Another source for the determination of long-distance matrix elements comes 
from fixed-target experiments. In Ref. |3l| Beneke and Rothstein find M^ 1S ^ = 30, 
having fixed < 0^ 2S \ 3 Si) > equal to 6.6 (in units of GeV -3 ). For the if)' resonance, 
the numerical values are M^ 25 -* = 5.2 having fixed < 0^ ls \ 3 Si) > equal to 4.6 (also 
in units of GeV -3 ) []. Therefore one gets R^ = 4.5 and R^ 2S ^ = 1.1, indicating 
again a smaller value for the 2S state, although subject to large uncertainties too 0. 

NRQCD matrix elements can also be determined from inclusive decays of B 
mesons. In a recent work [Q, Ma obtains a combination of < ( 1 5'o) > and 

< 0t {nS \ 3 Pj) > MEs in the framework of HQEFT, having fixed < of n5) ( 3 5i) > 
equal to 10.6(4.4)-10 -3 GeV 3 for the J/ip and ip' respectively. The resulting ra- 
tios (not exactly the same as in tables 1 and 2 since r takes the value of 1.13) 
are R^^ 1S ^ ~ R^( 2S ) ~ 2.3 which, although slightly smaller than those presented in 
the tables, are not inconsistent with them in view of the smaller r value and the 
uncertainties involved in the extraction [33|. 



7 Summary and last remarks 

In this paper we have carefully reviewed all steps to derive from the full QCD tree- 
level Lagrangian a complete transformed Lagrangian in terms of the heavy quark 
effective fields keeping the particle-antiparticle mixed pieces allowing for heavy 
quark-antiquark pair annihilation/creation. Let us note that such pieces are not 
generally neither used nor shown in similar developments in the literature, with a 



few exceptions |13, [14 



Indeed, it may seem quite striking that a low-energy effective theory could be 
appropriate to deal with hard processes such as QQ annihilation or creation. The 
keypoint is that, assuming a kinematic regime where heavy quarks/antiquarks are 
almost on-shell and moving with small relative momentum, the strong momentum 
dependence associated with the heavy quark masses can be removed, so that a de- 
scription based on the low frequency modes of the fields still makes sense, as in 
HQET applied to heavy-light hadrons. Such a kinematic regime can be well matched 
by heavy quarkonia and intermediate colored (QQ)s bound states predicted by the 
COM. 

In particular, we have focused on an annihilation process with initial-state quarks 
satisfying the Dirac equation of motion for free fermions. Thus, we have derived di- 
rectly from the &~ + > piece the heavy quark vector current coupling to a background 



7 In fixed-target experiments it is not possible to fit simultaneously both matrix elements because they display a 
similar shape as a function of the beam energy. 

8 In Ref. smaller values for M^ nS ^ are obtained as the inclusion of NLO corrections tend to decrease the values 
needed for the MEs. However, all values shown in tables 1 and 2 refer to LO extractions so we prefer not mixing them. 
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gluonic field, recovering a well known expression shown in the literature |20| allowing 
a non-relativistic expansion in the matching procedure of NRQCD and full QCD. 

In the second part of this work we have examined the meaning of using a HQEFT 
Lagrangian to account for heavy quarkonium production. Thereby we kept the £(~ + ^ 
and pieces in the effective Lagrangian, from which the £4 four-fermion piece 

can be "constructed" . The result looks formally the same as in conventional NRQCD 
although the fields are the HQEFT ones, i.e. giving self-consistency to all the 

framework. 

From a phenomenological point of view, we advocate the possibility of using the 
extended version of HQEFT for describing inclusive hadroproduction of charmonium 
states according to the color-octet mechanism. Indeed, the four different energy 
scales usually invoked in NRQCD (i.e. tuq, rnqv', rriQv' 2 and Aqcd) basically reduce 
to two in HQEFT: tuq and a typical hadronic scale, A, of the order of a few Aged's, 
in a reasonable accordance with the dynamics of charmonium systems. Conversely, 
inclusive production of bottomonium resonances should likely be better described 
by NRQCD because of the clearer difference among the scales. 

From inspection of tables 1 and 2, we can conclude that there are hints indicating 
that HQEFT may be considered as a candidate for describing J / ip and ip' hadropro- 
duction - although likely in a limited way because of the approximations in the scales 
and the lack of a static limit a la HQET when A/mq — > 0, as stressed elsewhere 
|5], Hence we give support to the hypothesis presented in Ref. of basing 
the power counting for charmonium hadroproduction along the lines of HQEFT, to 
some extent. 

As a final comment, let us remark that still the uncertainties in the phenomeno- 
logical extractions of the color-octet matrix elements (choice of numerical values for 
the heavy quark mass and renormalization/fragmentation scales, parton distribution 
functions, etc) are quite large. Moreover, new approaches have recently appeared 



(e.g. the krp factorization |34], |35]) while the experimental evidence is not yet conclu- 
sive to decide which framework should be the most adequate to describe hadropro- 
duction of heavy resonances. 
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